Zero Field Wigner crystal 
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A candidate for the insulating phase of the 2D electron gas, seen in high mobility 2D MOSFETS 
and heterojunctions, is a Wigner crystal pinned by the incipient disorder. With this in view, 
we study the effect of collective pinning on the physical properties of the crystal formed in zero 
external magnetic field. We use an elastic theory to describe to long wavelength modes of the 
crystal. The disorder is treated using the standard Gaussian variational method. We calculate 
various physical properties of the system with particular emphasis on their density dependence. 
We revisit the problem of compressibility in this system and present results for the compressibility 
obtained via effective capacitance measurements in experiments using bilayers. We present results 
for the dynamical conductivity, surface acoustic wave anomalies and the power radiated by the 
crystal through phonon emission at finite temperatures. 

PACS numbers: 



I. INTRODUCTION 



The question of the combined effects of disorder and 
interactions in correlated electron systems is one of the 
more important issues in solid state physics. The in- 
terest in this longstanding problem has been revived by 
the spectacular experimental results on the two dimen- 
sional electron gasi. From a theoretical standpoint, both 
disorder and interactions are challenges in themselves. 
By starting from the noninteracting limit, considerable 
progress was achieved in the understanding of the effects 
engendered by disorder in the system. The rationale be- 
hind such a limit stems from the expectation that not too 
strong interactions will lead to a Fermi liquid behavior, 
at least in high dimensions, and therefore, excitations of 
the system will behave as free fermions. For non inter- 
acting electrons in one and two dimensions, disorder is 
known to lead to the phenomenon of Anderson localiza- 
tion where, all states are localized 2 *^. Going from this 
non interacting limit to the interacting system has proven 
quite challenging. In three or higher dimensions, where 
Fermi liquid theory is expected to hold in the absence of 
disorder, a renormalization group studji^ indicated that 
disorder strongly enhanced the interactions hence push- 
ing the system away from the noninteracting point. In 
low dimensions, the situation is even more complicated 
because in this case, even in the absence of disorder, the 
interactions themselves have singular effects on the sys- 
tem, and can lead to non fermi liquid phases. For ex- 
ample, in one dimension, where the interacting system 
is known to be a Luttinger liquid^, the combined ef- 
fect of disorder and interactions have been shown to lead 
to behaviors quite different from that anticipated for a 
noninteracting systerrA 

Two dimensions is thus from this point of view doubly 
marginal: first, the effect of interactions even in the pure 



system is far from being completely settled and secondly, 
the disorder, although it leads to localization for the non- 
interacting system does it only marginally. It is therefore, 
natural to expect that interactions can alter this picture 
of marginally localized states. Consequently, the nature 
of the phase diagram of the two dimensional electron gas 
as a function of the disorder strength and interactions 
has been a subject of intense debate. There is however, 
one limit where the combined effect of interaction and 
disorder is amenable to an analytical study. This is the 
limit of strong interactions, where, the electrons in the 
pure system localize around the sites of a triangular lat- 
tice to form a Wigner crystaU2*ii. In this crystal, and 
far enough from melting, the particles become discernible 
by their position. Thus quantum effects are simple and 
manifest themselves in the quantification of the vibra- 
tions of the electronic crystal. When a disorder potential 
is added, its effect is to pin the elastic crystal^. The 
pinning of elastic structures, both classicaii&ilii^ and 
quantumiSiil, by disorder has been the subject of in- 
tense studies and various methods have been developed 
to tackle these problems. These tools have been used to 
study the magnetic field induced Wigner crystal forma- 
tion in the 2D electron gao ls i 19 i 20 i 21 i 22 . The optical con- 
ductivity and Hall coefficientsi2i2I were computed and 
found to be in excellent qualitative agreement with the 
experimental observations 2 ^ 2 ^. Combined with the the- 
oretical analysis, the optical conductivity provided clear 
evidence that for certain filling fractions, the phase real- 
ized in this system was indeed a Wigner crystal, weakly 
pinned on impurities. 

The formalism used to study the strong field crystal 
can be easily modified to obtain the properties of the 
B = crystal. This is of importance to experiments 
on ultra clean Ga-As heterojunctions and Si MOSFETS 
which exhibit the zero field metal-insulator transition. A 
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better way to parametrize the system is through the di- 
mensionless parameter r s — a/ as, where the lattice spac- 
ing a is related to the electron density n = 1 /ira 2 and clb 
is the Bohr radius of the electron in the 2D electron gas. 
One of the important questions is whether the insulating 
phase seen in these systems at low densities or large r s is 
indeed a disordered Wigner crystal. In the present paper 
we focus on analyzing the physical properties of a puta- 
tive Wigner crystal phase. Our approach complements 
the standard approach starting from the noninteracting 
limit. We investigate the effects of the disorder on such 
a system and compute various observable quantities. 

The plan of the paper is as follows. In Sec. [H] we first 
use a variational wave function to estimate the effective 
particle size i.e., size of the localized wave packet at a site. 
This parameter is expected to play an important role as 
will be detailed later in the paper. We then present the 
model and the variational solution. The remainder of the 
paper is devoted to the computation of various observable 
quantities: Sec. lIIII deals with the compressibility, Sec.HVl 
with the transport properties, in particular, the conduc- 
tivity, surface acoustic wave absorption and the power 
radiated by the pinned crystal, followed by a concluding 
section Sec.[V] Technical details have been relegated ton 
the Appendices. 



II. MODEL AND METHOD 

We will follow the same procedure outlined in Refill. 
I2TL In the crystal state, the wavefunction of the particles 
is localized around positions R^ . The wave function has a 
characteristic width £ centered around this site. Since the 
localized wavefunction renders the particles discernible 
by their position, the crystal can be considered as a col- 
lection of discernible particles of size £ that are labelled 
by their position R;. Then only the (quantized) vibra- 
tional modes corresponding to the displacements of these 
particles needs to be retained in any low energy descrip- 
tion of the system. Of course, this greatly simplifies the 
analysis and is the key ingredient which permits a res- 
olution of this problem. As pointed out previousl y 19 ! 21 , 
there arc three important lengthscales to describe the 
crystal: the lattice spacing a, the size £ of the particles, 
and the correlation length of the disorder rf which we will 
come back to later. For the case of a strong magnetic 
field the size of the particle is essentially the cyclotron 
radius 19,21 . However, for the B = Wigner crystal, de- 
termining the size of the particle £ is a question in itself 
since it is determined by the competition between kinetic 
energy and interactions. 



A. Wavefunction in the crystal state 

In this section, we use a variational wave function to 
determine this effective width £ as a function of the elec- 
tron density. We choose the ground state to be a Slater 



determinant of gaussian wave packets, whose width £ is 
chosen to be the variational parameter. The single par- 
ticle gaussian wave packets centered around a site i with 

coordinate Ri, take the form ipi(r) — wfexp- 



(r-Rj 



In the absence of a magnetic field, the total variational 
energy per site is given by 



E(0 



m£ 2 2e£ 



(1) 



where the first term is the kinetic energy of the electrons, 
and the second term is the potential energy arising from 
coulomb repulsion and 



V(Ri) 



R 2 



R 2 



Io is a modified Bessel function of the first kind and R; 
denotes the lattice sites, m, e and e denote the mass, 
charge and dielectric constant respectively. Minimizing 
the energy E with respect to £ results in a self-consistent 
equation for £. Excepting the limit of low densities, where 
one can obtain an analytic expression for £, this varia- 
tional equation is rather complicated to solve analytically 
for arbitrary densities, thereby requiring the use of nu- 
merical techniques. In this paper, we discuss only the low 
density limit or equivalently the deep crystalline phase. 
In this limit, as the density is lowered, the effectiv e dis- 
tance between neighboring sites increases as a = y/l/im 
and the effective width of the wave packet increases in 
manner such that a/£ 3> 1. Therefore, it is sufficient to 
retain only the term Rf = a 2 in the expression for V . 
Expanding the exponential in £/a, we obtain following 
expression for the energy for low densities: 



E 



m£ 2 



C 



e 2 e 



(3) 



The constant C depends on the coordination number of 
the lattice and the number of terms retained in the sum. 
Retaining only nearest neighbor sites on a triangular lat- 
tice, we obtain C ~ 1.88. Minimization of E with respect 
to £, yields the result, 



,-3/8 



(4) 
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an = is the Bohr radius of the electron. Note that 

' me 

despite the fact that £ increases with decreasing density, 
the ratio £/a ~ n» decreases with decreasing density in 
accord with the assumptions made above, thereby justi- 
fying the expansion © . We emphasize that this result is 
valid only deep in the crystalline phase. Close to melting, 
exchange terms become important and we need to solve 
the equation with the full V to obtain a reasonable varia- 
tion of the effective particle size with density. Moreover, 
when exchange terms become non-negligible, the very hy- 
pothesis of having discernable particles collapses, and an 
effective theory different from the one presented in the 
ensuing section is required. 



3 



In the presence of a magnetic field, the same calcula- 
tion can be done to obtain £ as a function of the density 
and the field. For the purely gaussian wave packet, the 
magnetic field contributes a term mo; 2 £ 2 /2 to the en- 
ergy per site. In the limit of ultra strong magnetic fields 
dominating the coulomb repulsion, we recover the result 

that £ = which is just the cyclotron length. For 

arbitrary fields, 



m 2 u 2 + Cme 2 /ea?\ 

However, the very assumption of a gaussian wave packet 
implies that the electrons are confined to the lowest Lan- 
dau level and hence the above result holds, strictly speak- 
ing, only for strong magnetic fields, and should be viewed 



as a convenient interpolation formula. 



B. Elastic Hamiltonian 

Now that we have the important parameters charac- 
terizing the crystal, we present the elastic Hamiltonian 
describing the crystal phase. We use the same recipe 
as outlined in Ref. Il9ll2ll and recall only the main steps 
here. In the crystalline phase, the electrons occupy the 
sites of a triangular lattice with a lattice constant a which 
is related to the density of electrons by n ~ (7ra 2 ) -1 . A 
particle at a site i is displaced from its mean equilibrium 
position denoted by Ri, by u(Rj,i). In the continuum 
limit, the vibration modes of the crystal lead to the fol- 
lowing elastic action in Fourier space 



S = X! J d2q [^(q, ^n){Pm^n + £V(q))u T (-q, -w n ) + u L (q, io n ){p m u? n + fi L (q))u L (-q, -u n )] 

+ I d 2 rdtV(r)p(r) (6) 



where the transverse (longitudinal) displacements ut 
(ul) are related to the cartesian displacement u as fol- 
lows: 

u a (q) = w L (q)q Q + u T (q)e a/ 3q / 3 (7) 

q = q/|q| is the unit vector along q and e a p is the fully 
antisymmetric tensor with e xy = 1. These two modes 
can be interpreted as the shear (ut) and the compression 
(ul ) modes respectively. The Matsubara frequencies are 
u) n = with (3 = 1/T being the inverse temperature. 
The terms quadratic in the frequency represent the ki- 
netic energy of these modes. p m — -^,p c = -^z are 
the mass and charge densities respectively. The elastic 
energies of these modes are giverj 2 ^ by 

^i(q) =di|q|+c L q 2 

r> f \ 2 ^> 

f2 T (q) = c T q 

where cl,ct,cL are elastic constants and, eo is the di- 
electric constant of the substrate. For the classical crys- 

2 

tal on the triangular lattice, one has 2 ^ cl = —0.18-^, 

2 2 

ct = 0.04— and d — — . These values for the elas- 
tic modulii in JHJ, are valid only for low densities, deep 
in the crystal phase. For arbitrary densities, the elastic 
constants depend on the scale £ and change drastically as 
one approaches melting. This point should be taken into 
account while comparing theoretical results for the den- 
sity dependence of various quantities, with experiments. 
Unfortunately a rigorous estimate of the elastic constants 



as a function of the density is still lacking. The linear q 
dependence in the compression mode fii(q) arises from 
the coulomb repulsion. This is because a longitudinal 
deformation changes the density profile of the crystal 
which costs Coulomb energy. On the contrary, trans- 
verse modes can be excited without changing the density 
and are thus purely elastic. 

Finally, the last term in the action describes the cou- 
pling to disorder, modelled here by a random potential 
V. The density of the particles p(r) is given by 

p(r)=^5(r-R i -u i ) (9) 

i 

where 8 is a 5-like function of range r/. Since the disorder 
can vary at a lengthscalc rf which can a priori be shorter 
or comparable to the lattice spacing a, the continuum 
limit u(Rj) — > u(r), valid in the elastic limit |u(Rj) — 
u(Rj + i) <ti a, should be taken with care in the disorder 
termSSiSi. This is done using the decomposition of the 
density in terms of its Fourier components 

p(r) ~ po - PQ V ■ u + po e iK < r -»W (10) 

where po is the average density and K are the reciprocal 
lattice vectors of the crystal. 

The above formula shows that one should distinguish 
between various parts of the disorder. Disorder that 
varies at lengthscales much larger than the lattice spacing 
couples only to the /9qV • u term. Although such a term 
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can modify the structure of the crystal quite strongly, it 
does not lead to pinning and disappears from the trans- 
port properties. Pinning, and thus the dominant contri- 
bution to transport coefficients, is however achieved by 
the disorder that varies on lengthscalcs comparable to 
the lattice spacing. In hctcrojunctions, such a disorder 
is expected to be present due to interface roughness aris- 
ing from epitaxial growth of the semiconducting layers 
sandwiching the electron gas. The real system is further 
complicated by the presence of a long range disorder po- 
tential arising from the presence of ineffectively screened 
dopants outside the plane of the electron gas. However 
such a disorder varies slowly compared to the scale of the 
lattice spacing, and we will thus assume that we can ne- 
glect it as far as the transport properties are concerned. 
Moreover, if the disorder is weak and leads to collective 
pinning one can assume a gaussian distribution for V, 
since each volume of the system will average over a large 



number of independent impurities. In this case, 



V(r)F(r') = A r/ (r-r') (11) 
where A r/ is effectively a delta function with a range r/. 



C. Gaussian variational method 

We treat ®) using a variational methodi^^. We 
present here only the main steps of the treatment and 
refer the reader to Ref. |^] for details. Many of the tech- 
nical details and subtleties of the m ethod can be found in 
the literature (see e.g. Ref. IH 123 for a review) . We first 
average (|SJ) over disorder by introducing replicas. This 
averaging results in an effective action which involves in- 
teractions between the n replicas, given by 



s =1^2 f u t(^ u n )[p m ujl + n T (q)]u£(-q, -w n ) + <(q, u n )[p m wl + fii(q)]<(-q, -ui n ) 

- ^ [ d 2 r f f drdr' V A K cos [K • (u a (r, r) - u\v, /))] (12) 

1 J J ° J ° a,b,K 



The replica indices, a, b run from 1 to n. The size of 
the particles £ and the finite correlation length of the 
disorder restrict the sum over K to values of K smaller 
than -Kmax ~ t/ max(r^ , £). Ak is taken as a constant 
Ak = A for K < K max and zero otherwise. The physical 
disorder averages are recovered in the limit n — > 0. 

We now search for a variational solution to (|12fl by 
using the best quadratic action approximating 112(1. We 
use the trial action 

S ° = YB / £ <(q^n)(G- 1 )^(q,^ l )<(-q,-^) 

(13) 

where the whole Green's function (G 1 )°*(q, w n ) are 
variational parameters. The variational free energy is 
now given by 



F-vax = Fq + (S — So) So 



(14) 



The pertinent solution of these saddle point equations 
(|15f) solved in the limit of the number of replicas n — > 0, 
breaks replica symmetryiiiSi* 2 ^. 

The connected Green's functions defined as (G, 
X^(G -1 )^ for the two modes are given by 



c l^v 



G cT {q, iu) n ) 
G cL (q, iu n ) 



1 



p m tol + f2r(q) + I(w n ) + £(1 - S n ,o) 
1 

p m u;l + fii(q) + I(iu} n ) + E(l - 5 nfi ) 

(16) 



where the disorder induced pseudogap is found to be 



E = c T R- 2 = c T (27r 2 )-6i?- 2 (a/^o) 6 



(17) 



The variational parameters are then determined by the 
saddle point equations 



0(G-i)«*(q,w„) 







(15) 



2 __ 

with £ = max [ r /)^]j c t — ^-f— and R a — ct/to 2 vA. 
As explained in Refs. l2lll23 . the function I(iu> n ) is deter- 
mined by the following self-consistent equation derived 
in the semi-classical limit: 



5 



1 



ft T (q) + £ «i(q) + S 



fir(q) + p m w^ + /(iw n ) + £ Qz,(q) + + J(iw n ) + s . 



(18) 



In the limit of small Lu n , the above equation can be solved 
analytically to obtain 



'2p m S[l 



4c T S , , <P 



(19) 



This equation for / can be continued to real frequencies 
and then solved numerically. The solution of the vari- 
ational equations as obtained above thus allows us to 
extract numerous physical quantities of the disordered 
system. We will first discuss the compressibility in the 
forthcoming section and the transport will be discussed 
in Sec. |IV| 



III. COMPRESSIBILITY 



is easily computed to be 



k oc lim — , 

q ^o,uj n =o ojI + Sl L (q) 



(21) 



Since ^(q) oc q, the compressibility is zero. This arises 
from the fact that the change of density induced by p 
occurs for a constant neutralizing background. The com- 
bined system of the crystal plus the background thus be- 
comes charged which costs too much energy and inhibits 
this kind of density fluctuation. 



Pi 



-nil 



In this section we discuss the compressibility of the 
pinned Wigner crystal. Although the compressibility 
is usually simply related to density-density correlations, 
here one should take special care in defining this quantity 
due to two complications inherent to the Wigner crystal: 
(i) for charged systems it is crucial to know whether one 
keeps the system neutral or not when letting the den- 
sity fluctuate, (ii) since in presence of disorder one ex- 
pects glassy properties, the question of the timescale over 
which the measurements are performed is relevant. 

The naive way to define the compressibility K would 
be to add to the Hamiltonian the standard chemical po- 
tential term 

H = J d 2 r{p{v) - p ) (20) 

and to compute the average change in density (p(r) — 
pa). In linear response, the compressibility is given by 
the density density correlator. Using (J7J and (|10H the 
compressibility is directly related to the correlation of 
the longitudinal displacements. For the pure system this 

I 



FIG. 1: Capacitance measurement, which gives access to the 
compressibility of the system. A voltage difference /i is ap- 
plied to a capacitor. Here for simplicity, the capacitor is made 
of two planes of the 2DEG. 



To have a non-zero compressibility, even for the pure 
system, it is thus important to change the density of 
the Wigner crystal, while at the same time changing the 
background to maintain neutrality. An experimental way 
of doing this is to perform a capacitance measurement 30 . 
We show here how to compute such a capacitance for 
the case of two Wigner crystals planes^i separated by a 
distance d, as shown in Fig. ^ Though the experimen- 
tal setup (two layers and a gate) considered in Ref. is 
slightly different from the one considered here, we expect 
the results obtained here to hold for this experimental 
geometry as well, since the crucial ingredient, namely 
ensuring the neutrality of the system is preserved. 

The Hamiltonian of the system depicted in Fig. is 
thus 



H = H° 1+ H° + - £ / V aP (r-r')[p a (r)-p ][p p (r')-p } + L / [ Pl (r) - p 2 (r)] (22) 

(a,/3) = l,2^ r ' r ' Jr 

I 

where H± 2 are the elastic Hamiltonians of each 2DEG excluding the Coulomb interaction. V a p is the usual 
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Coulomb interaction in and between the planes. If one 
assumes that the system is neutral in the absence of /i, 
then within linear response, the charge on one plane when 
a potential p is applied is 



<Pl) = |[(PlPl> ~ (PlP2>] 



(23) 



Thus lj2*3)l gives directly the capacitance C = 
More generally, we define the capacitance for finite q and 
uj (e.g. corresponding to a time and space dependent 
potential p) 



CV,q) = -[(pipi)u,, q - 



(P1P2) 



(24) 



where the subscripts u>, q indicate that the correlation 
functions are evaluated for finite momentum and fre- 
quency. For a metal, a RPA evaluation of (|21f) gives 
back^i the standard formula for the compressibility of a 
charged system. We recall the calculation in Appendix 1X1 
for convenience. One can use this general formula to com- 
pute the capacitance for the Wigner crystal. Since one is 
interested in the limit of small q and u> one can use the 
decomposition of the density JTDJ. The Hamiltonian © 
becomes 



H = H 



m 



P 2 o 



2- £ 

(«,/9)=l,2 



^(r-rO[V.u Q (r)][V- U/3 (r')] (25) 



where Hf r 2 is the part of the Hamiltonian that does not 
contain the long range (i.e. for q ~ 0)) of the Coulomb 
interaction. Ui and U2 denote the displacement vectors 
in the two planes. Since Vii(q) = V2 2 (q) oc l/|q|, the 
third term in l|25|l obviously yields the part proportional 
to q in the bulk modulus for an isolated plane (see JHjl). 



Note that in writing l|25|) we have neglected the short 
range part of the interplane interaction 



Hti = I V a0 {r-v')pl J2 ^ 

Jr ' r ' K,K'#0 



(r-ui(r)) e -iK'-(r'-u 2 (r')) 



The translational invariance restricts the sum to K = K' 
terms. This essentially contains the Fourier transform 
of the interplane potential which using (pTTJ) is found to 
behave as 



V 12 (r)e 



iK-T 



(2ir)e 



-Kd 



K 



(27) 



Thus if the planes are at a distance d much larger than 
the lattice spacing a of the WC, this term is obviously 
much smaller than the long range part of the intraplane 
interaction, and can be neglected. It can also be ne- 
glected if one of the planes is an homogeneous electron 
gas (i.e. a simple metal). In the case where the two 
planes are close enough this term should be retained and 
can lead to interesting effects such as the locking of the 
two Wigner crystals together 3 ^. The in-plane coupling of 
the higher harmonics of the density contribute to H SR 
and generate the regular non-singular part of the hamil- 
tonian (i.e. the part proportional to q 2 in the elastic 
coefficients). In fact thi s ap proach is synonymous with 
the method used in Ref. 25 to calculate the elastic coef- 
ficients. The derivation is done in Appendix iBl 

Let us now use the general formula l|24|l to compute 
the capacitance of the pure system. From (fltjfl , the long 
wavelength part of the density is given by 



Mq) 



-ipoqu aL (q) 



(28) 



One thus needs only the longitudinal part of the action 
to calculate the compressibility. 



J 



S : 



<(q) 



,(q) 



^ s L r (q) + p^n(q) 



r 



pll 2 V l2 {ai) 
fi B L r (q)+p§<7 2 Vn(q) 



4(-q) 
d(-q) 



(29) 



where S^ r (q) is the "short range" part of the elastic co- 
efficients |jSJ. Using (|2*51l and (|2*4l one obtains the capac- 
itance 



C(wj n ,q) = ~ 



1 2 

Pol 



2 p m uj\ + nf(q) + plq^M) - Vi 2 (q)] 

(30) 

The thermodynamic compressibility is given by the value 
of C for uj n = 0. Note that in this case one obtains the 
same value by considering the retarded correlation func- 
tion (doing the analytic continuation iu n —* to + iS) and 
taking the limit u — > first (for a fixed q). The diver- 
gence arising from the long range part of the Coulomb 



potential now cancels since 



Vn(q) - Vi 2 (q) = / d 2 re^ r [- - -=L 



Vr 2 + d 2 
(2?r)(l - e-i d ) 



(31) 



is non divergent when q — > 0. This of course traduces 
the fact that the global system remains neutral when the 
potential is applied. Using (|31|l one obtains 



C{u n = 0, q -> 0) 



2 2 

PoT 



2nf(q) + 4nd(p 2 q^ 



(32) 
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One thus recovers that the inverse capacitance is the sum 
of a purely geometric term and an electronic one (see 
(I A4|) ) . The electronic capacitance is simply given by 



C e 



lim 



pfa 2 



pI 



'(g) 2 CL 



7TT- « "37 (33) 



This term, which for a normal metal is simply related to 
the static compressibility, is thus negative for a Wigner 
crystal, if one uses the classical estimates 2 — for the elastic 
coefficients. The fact that a system of discrete charges 
can lead to such effects has been noted before for classical 
Wigner crystals (see e.g. Ref. |33| and references therein) . 

Let us now turn to the disordered case. In the pres- 
ence of disorder, the same calculation can be repeated 
to obtain the compressibility. To do this, we assume 
that the disorder potentials in each layer is drawn from 
independent distributions and repeat the variational cal- 



culation of the previous section. The trial action H13|) 
has now two additional indices that denote the two lay- 
ers. For simplicity, we will only denote these indices here, 
all the others being implicit. The trial action is thus (in 
these indices) a two by two matrix with the Green's func- 
tions G~l where p,, v = 1,2. Obviously G\\ = G22 and 
G12 = G21. From Ij24(l and the inversion of the two by 
two matrix the capacitance is thus given by 



C7(w,q) 



pW 



[G11 — G12] 



pfo\ 



1 



2 l G^-Gl 



(34) 



where all other indices (replica, longitudinal, q, lu) are 
implicit. The variational procedure can now be repeated 
to determine G\\ and G12. Because the disorder is in- 
dependent from plane to plane, the trial free energy has 
the following matrix structure 



J 



\p m u> 2 n + nf(q) + Poyyii(q)]Gn + F[Gj 
p 2 q 2 V 1 2(q)G 2 i 



P§9 2 V r i 2 (q)Gi2 \ . 

[p m oj 2 n + nf((q)+pl, q 2 Vn(ci)}G22+F[G 22 \ ' 

I 



where F is the same as the term induced by disorder 
averaging in a single plane. The minimization l|15(l now 
gives the self consistent equations 



G n X = [p m U 2 n + nf(q) + P6V ^n(q)] + F'[Gn] 



(36) 



Gu 1 = P%q 2 Vi 2 (q) 
Note that the in-plane inverse Green's function is iden- 



tical to the one in the absence of interplane interactions, 
and thus is given by the variational solution of the pre- 
vious section. The interplane inverse Green's function is 
trivial and simply given by the long range part of the 
interaction between the planes. 

In the presence of disorder, the compressibility is given 
by the connected Green's function! 6 . Using (|34|) . 136(1 
and <|lf>|> one obtains 



G(w„,q) 



PW 



2 p m ujl + n s l(q) + P6V fMq) - V 12 (q)] + + E(l - S„,o) 



(37) 



r 



As for the pure system, we see that the geometric capac- 
itance contributes to the total capacitance. Isolating the 
electronic part of the capacitance we find 



G ei (w„,q) 



pW/i 



ft s /(q) + /^„) + £(l-<5„,o) 



(38) 



(|38|l thus shows that in the presence of disorder one has to 
distinguish between the thermodynamic capacitance and 
the dynamical one. If one considers the thermodynamic 
capacitance, then should be computed for u n = 0. 
In that case all contributions from the disorder disappear 
in (|38|l and the capacitance is identical, within the varia- 
tional approximation, to that of the pure system. On the 
other hand, if one considers a capacitance measurement 



in response to a modulation at finite frequency, which 
is certainly the case experimentally, one has to take the 
analytic continuation iuj n —tuj + iS first. In this case one 
obtains 

PoV/2 



G e ;(w,q) 



(39) 



Q s /(q)+7(a;)+E 

As for the pure case, one should take the limit u> — > 
first. In this limit, the function I(u>) is regular and tends 
to zero. Taking q — > 0, we see that since the mass term 
S prevails in the denominator the compressibility is zero. 
This is a consequence of the fact that the pinning by 
disorder renders the system inflexible to charge modu- 
lations. Note that contrary to the pure case, the ther- 
modynamic response and the slow dynamic one are not 
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equivalent in the disordered case. This is not surprising 
considering the glassy nature of the system. 

As for the effect of a magnetic field on the compress- 
ibility, naively one would expect that it has no effect since 
the compressibility is a static property. This is what is 
found within the variational approximation. However it 
is important to note that, although there is no explicit 
magnetic field dependence of the compressibility, there 
is still a variation of the thermodynamic compressibility 
with the magnetic field accruing from the field depen- 
dence of the clastic constants of the crystal. 



Rea o.4 




n=1 

n=0.9 

n=0.8 

n=0.6 

n=0.4 



IV. TRANSPORT 

We now focus on the consequences of disorder for the 
zero temperature transport properties. This is of partic- 
ular importance to the materials which exhibit the zero 
field metal-insulator transition where transport is indeed 
the main probe of the physics. Up to now, the experimen- 
tal emphasis has been on finite temperature resistivity 
and magneto-resistance measurements. However for the 
Wigner crystal such quantities are difficult to compute 
theoretically since they are dominated by the defects in 
the systemSi. It is thus difficult to use them as a probe 
of the Wigner crystal nature of the underlying phase. As 
explained in Ref. Il7ll2ll the optical conductivity does not 
suffer from such a problem and can be reliably computed 
theoretically from the elastic Hamiltonian of the crystal. 
Here, we present results for the dynamical conductiv- 
ity, surface acoustic wave measurements and the power 
radiated by the crystal. These measurements done in 
the quantum Hall samples at high fields were instrumen- 
tal in clarifying the physics of these systemaS&Sii&SiSi. 
Clearly similar measurements are called for in the zero 
field samples, and would be crucial to understand the 
physics of the insulating phase. 



A. Conductivity 

The conductivity of the disordered crystal can be ob- 
tained from the displacement-displacement correlation 



functional and is given by: 



cr a p{u) = ip 2 c uiG al3 {q = 0, co + ie) 



(40) 



where (iv = x,y and G MIy (q,w) = (w^(q, co)u u (q, co)) 
are the displacement Green's function. Since, there is 
no magnetic field, the longitudinal resistivities a xx (ui) — 
a vv {co) = <j{u>). These functions are related to the con- 
nected Green's functions of l|16|l and using an analytic 
continuation ico n — > co + ie, we obtain 



S + I(co) 



(41) 



As seen in other disordered elastic systematise, the con- 
ductivity is completely determined by £ and I {to). For 



FIG. 2: The real part of the conductivity <j(lo) (in units of 
e 2 1 ' yfA^in ) as a function of the frequency measured in units 



of \ 2l for the regime where the disorder correlation length 

Tf < £ for various electronic densities n(in units of (-Ka%)~ 1 ). 
These curves have been obtained using the classical values of 
the elastic moduli and the parameter A2 is given in the text. 

the pure crystal, £ = I(to n ) = 0, and one recovers that 
the real part of the conductivity exhibits a co — Drude 
peak and zero finite frequency conductivity, while the 
imaginary part varies as 1/co. In the presence of disor- 
der, the crystal is pinned collectively forbidding any slid- 
ing motion of the crystal. Consequently, the Drude peak 
is annihilated and the dc conductivity Kea(uj = 0) = 0. 
This results in a non-zero finite frequency conductivity 

and the appearance of a peak at a new scale co p ~ \[^F~ 

called the pinning frequencji^i. The pinning frequency is 
thus determined by l|17jl. The pinning peak is broadened 
by the dissipative term I(co) which also generates, as can 
be seen from 1)41 [I. a slight shift of the pinning peak from 
lo p due to the real part of I(co). Using l|19[) . we find that 
for co <C co p . 



Recr(w) ~ p 
Imcr(u') 



4c T S 



(P 



-'2p m E[l + ^log(-^)]^ 



, CO 

V CT S' J S" 2 



(42) 



while for high frequencies co 3> co p 



Re<r(a;) 



Pc 



Prn^ 3 



(43) 



The behavior of the conductivity for all intermediate val- 
ues of co can be obtained by solving numerically l|18|) and 
using the solution for I(co) in (|41(l . Referring to (|17|l . we 
note that the result for T, depends on whether the disor- 
der correlation length is larger or smaller than the size of 
the particles of the crystal (r/ £) . This has serious im- 
plications for the density dependence of the conductivity 
as can be inferred from Figs. [21 and [3] where, we plot the 
conductivity for different values of the density in the two 
regimes rf < £ and rf > £. 

The density dependence of the pinning frequency and 
the height of the pinning peak can now be evaluated in 
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FIG. 3: The real part of the conductivity <j(uj) (in units of 
(? I \J A\m ) as a function of the frequency measured in units 

of Jj± for the regime where the disorder correlation length 

Tf > £ for various densities n ( in units of (■kcl 2 b )~ 1 ). As 
before, these curves have been obtained using the classical 
values of the elastic moduli and the parameter A\ is given in 
the text. 

a straightforward manner. The peak height is given by 
the expression 



£ i 
P = <j(u} p ) = p*J—Im- 

\ Pm. A 



(44) 



If we neglect corrections to the pinning frequency arising 
from /, and assume, then I(to p ) = T,Z where Z is some 
complex number, and in this case, the leading density 
dependence of the peak height is given by 



P oc 



(45) 



The behavior of the peak height is thus entirely deter- 
mined by £. To obtain the actual density dependence, 
one needs to know the elastic moduli of the crystal. Us- 
ing the results of Ref. l25l we find that in the regime, 
r/ > £, the electrons see the bare disorder and 



£ = Aya 



Rc 



be 2 



where A\ = (27r 2 )"5 Ae/by/ire 2 r 6 f and b 
classical triangular lattice. This leads to 

_ 3 

LOp oc n 4 
P oc 



(46) 



0.04 for the 



(47) 



Consequently, the peak position shifts to higher frequen- 
cies and the peak height decreases with decreasing den- 
sity as is seen in Fig. |3] 

In the opposite regime where the disorder correlation 
length rj < £, the effective particle size being larger, the 
particle sees an averaged effective disorder. Using the 



classical values for the elastic moduli, we find that 
£ = A 2 n^ 



Rr.= 



be 2 



(48) 



A-zey/n 



where A 2 = (2n 2 )-i keirl (C/a B )^/be 2 and C = 1.88 for 
the triangular lattice cf. |3| Hence 



P cx n» 



(49) 



This implies that as density is decreased, the pinning 
peak shifts to lower frequencies accompanied by a con- 
comitant decrease of peak height as shown in Fig. [21 The 
density dependence of the pinning peak is a test that can 
prove whether the insulating phase is indeed a Wigner 
crystal. 

Another interesting quantity is the threshold electric 
field for depinning of the crystal. This field is related^ 
to the parameters of the problem by the relation Et — 
ctRc 2 £, where £ = max[r/,£]. Using ^46148(1 . we obtain 



Et oc n 
Et <x n~ 



r f >£ 



(50) 



Note that for r/ > £, the threshold field increases with 
inverse density, implying that the crystal gets more and 
more pinned in the low density regime. However, since 
f increases when the density decreases according to (@J 
on has to crossover to the other regime r/ < £. Depend- 
ing on the correlation length of the disorder, the system 
can thus cross over from one kind of dynamical behav- 
ior to another as the density is varied. In the first case, 
the pinning frequency and the threshold field increase 
with decreasing density. When the particle size exceeds 
the length scale of disorder, the pinning frequency and 
the threshold decrease with decreasing density. In both 
cases, the peak height decreases with decreasing density. 
This is compatible with the reduction of spectral weight 
with decreasing density. Such a crossover would lead to 
a maximum in the threshold field as the density is de- 
creased. 

If we now turn on a magnetic field perpendicular to 
the plane of the crystal, the result for E is unchanged 
as long as the magnetic length l c > £. For weak fields, 
the peaks in the diagonal conductivity now occur at the 

2 

frequencies aip 3 = w c ± ^p =0 ± 8 Jb=o ■ F° r strong enough 
fields such that l c < £, we revert back to the usual strong 
field physics discussed in Ref. l2ll 



B. Surface Acoustic Wave Measurements 

Another physical property that can be extracted from 
our results pertains to surface acoustic waves (SAW). 
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SAW measurements were used in the past to obtain the 
crystal dispersion relations^. In SAW measurements, 
the system is excited by a finite q SAW with a frequency 
uj = vq/2ir, where v is the velocity of the SAW. As the 
excitation traverses the system, the SAW propagation 
is completely affected by the piezoelectric interaction re- 
sulting in an attenuation and a shift in the SAW velocity. 
This shift dv and the attenuation rate n are essentially 
determined by the density response function of the sys- 
tem and are given by^i 

Av _ a 2 1 

V 2 1 i [£Ml2 

i + [ --J ( 51 ) 
_ qa 2 <j{g)a m 
K 2 al+a(q) 2 

where a is the piezoelectric coupling constant which de- 
termines the interaction between the SAW and the elec- 
tron system and a m = v(e + eq) where e and eg are the 



dielectric constants of the medium and vacuum respec- 
tively, and the momentum dependent longitudinal con- 
ductivity of the crystal a(q) = Rea(q,uj — vq). Notice 
that Av and n provide an indirect measurement of the 
conductivity at finite q. Using the results of the preced- 
ing sections, we obtain 

a(q, lo) = % G L (q, oj) + ^-G T (q, oj) 52 

q A 1 

where 

Gl(<Z ' W)= X + n L (q)-p m OJ 2 + I(oj) 

I ( 53 ) 

Gt(<7 ' w) = s + n r (g)-Pm^ 2 + /H 

Using I(uj) = h(u)) + il^u) one gets 



Recr(<7, u — vq) 



QyHvq) 



[E + q 2 (c L - p m v 2 ) + dq + h(vq)} 2 + h(vq) 2 [E + q 2 (c T - p m v 2 ) + h(vq)} 2 + h{vqf 



where I\ and I2 are the real and imaginary parts of I{oS). 
For small frequencies since IzioS) cx u>, and S is the dom- 
inant scale, 



Recr(<2, lo = vq) = T 



v 2 q 2 



(55) 



where T = ^/2p m E[l + ^log(^)] is a density de- 
pendent parameter independent of q. Consequently, one 
has 



Aw 

v 



1 



1 



rlVV 



l£ 2 o 

Tv 2 q 3 j: 2 a. 



(56) 



2 E 4 ct^ + T 2 v 4 q 4 



Using our results for I(u>) and the wave vector dependent 
conductivity, we can thus calculate the shift and attenu- 
ation for all values of q. The results are shown in Fig. 21 
In both the regimes r/ 5* £i despite a reduction in the the 
magnitude of the SAW anomaly with decreasing density, 
the width of anomaly is expected to be strongly density 
dependent and is essentially dictated by the competition 
of the disorder scale r/ and the effective particle size. 
These results provide an impetus for detailed SAW ex- 
periments in the 2d electron gas systems. We reiterate 
that the above derivation has been done with the classi- 
cal values of the elastic constants, and is thus valid deep 
in the crystal phase. A full quantitative prediction would 



Av/v 0.8 




FIG. 4: The typical profile of the SAW velocity shift (in arbi- 
trary units calculated for a m = 1 and o? = 2 ) as a function 
of the frequency to — vq and density. In both the regimes 
of disorder, the position of the minimum shifts to higher fre- 
quencies and the depth decreases as density is decreased. The 
quantitative behavior will be governed by the value of E. cor- 
relation length tf > £. The curves are plotted for various 
values of the density n. 



require a precise knowledge of the variation of the elastic 
constants with density, a piece of information missing at 
the moment. 
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C. Power radiated by the 2d crystal 

Another possible probe of the physics of the electron 
gas concerns thermometry. This involves a study of 
transfer of heat generated by a current in the electron 
gas to its three dimensional environment. This occurs 
via phonon emission at any given temperature and is re- 
lated to the internal dynamics of the electron gas. These 
measurements were done on the electron gas both in the 
absence and presence of magnetic fields 40 - 41 ' 42 . The ex- 
pression for the power radiated by the 2d gas at a tem- 
perature T into cooler environment, considered to be at 
zero temperature is given by^i: 

P = ^c(Q)|Af(Q)| 2 n s (/3^(Q))G(g, W (Q)) (57) 
Q 

where G(q,ui) cx wRe(l/cr(<7,a;)), q is the projection of 
the phonon momentum Q onto the plane of the WC, M 
is the electron phonon matrix element and nn is the usual 
Bose distribution function. Following Ref . l4fi we assume 
that the electron phonon coupling is dominated by the 
piezoelectric effect and | M(Q) | 2 ~ 1/Q. In the absence of 
any external field and for high mobility systems, G cx uj 2 
for ui = w(Q) oc Q and and the predicted power radiated 
P cx T 5 was found to be in accord with experiments. 
Based on power counting arguments, it was also shown 
that the radiated power in the IQHE phase F a T 4 
In the light of these experiments, an obvious question is 
what is the power radiated in the pinned crystal. A naive 
approach would be to suppose that finite temperature 
contributions to Rel/cx are small at low temperatures, 
hence permitting one to approximate G(q,u>) by its zero 
temperature value. In this case, we find that for the 
frequency range < to < uj Pi Rel/crfq, uj( q)) doe s not 
vary significantly. Following the lines of Ref. I40l4ll since 
the typical values of wave vectors and frequencies of the 
phonons are small compared to the scales over which the 
inverse conductivity changes, we replace G (x uj. Again a 
simple power counting results in P ~ T 4 akin to what is 
seen in the IQHE phases. 

A real calculation of P requires a knowledge of the 
finite temperature conductivity. At low enough tem- 
peratures, topological defects are thermally excited and 
they constitute the primary contribution to the ultra- 
low frequency conductivity. Though our formalism can- 
not treat the defects, we present heuristic arguments to 
justify our above result. A reasonable hypothesis for 
the low temperature frequency dependent conductivity, 
is a(q,u>,T) ~ a(q,uj,0) + ad(q, w,T), where ad is the de- 
fect contribution. The presence of ad leads to a finite but 
small dc conductivity a(T) in the long wavelength limit. 
We expect this to be the dominant contribution for small 
frequencies to < LUd and for frequencies u) > u>d, we ex- 
pect the finite frequency contribution to dominate. The 
cut-off frequency LUd is fixed by a(u>,T = 0) ~ ad(T). 
Since for small frequencies, a(uj,T = 0) ~ Aui 2 , we 
obtain u>d — \J &d{T) /A. This leads to the form that 



Re(l/a(q,u,T) 



7m 



Re(l/a(q,u,T = 0). Using 



this form in (|57|l . and doing the usual power counting, 
we obtain P cx T 4 + aT \J ad(T). Since in the limit 
of low temperatures <7<j(T) is expected to be exponen- 
tially small, the power spectrum is dominated by the T 4 
term, thereby justifying our naive analysis of the preced- 
ing paragraph. This coincidence is related to the fact 
that both the pinned crystal and the IQHE phases are 
incompressible. Clearly, the variations seen in the in- 
verse conductivity will modify this result, but we expect 
the result to hold for some range of low enough tempera- 
tures with a crossover to a different scaling form at higher 
temperatures. 



V. CONCLUSIONS AND EXPERIMENTAL 
PROPOSITIONS 

Before we conclude, we delve into the possibility of ob- 
taining quantitative estimates for the theoretical scales 
R a and E in the so called WC phase. These are very 
important for a real comparison of theoretical predic- 
tions with experimental data. The available experimen- 
tal data which includes mobility and resistivity measure- 
ments, permits one to calculate basic quantities like the 
the fermi momentum kp and scattering rates. The fermi 
momentum can be expressed as a simple function of den- 
sity 7ra 2 = n or r s = a/as, the mass renormalization 
factor m r and the relative permittivity e r of the mate- 
rial. The Bohr radius of the electron in the given material 
is as = e r /m r Q.53A° . Using the relation for the Fermi 
energy of the electron gas, Ep = , 2 = h 2 k 2 F /m* , one 
finds 



1.768n5 



■1. 



x 10 



10 



e r r s 



(58) 



kp is expressed in units of inverse meters. We also find 
that the plasma frequency of the electron gas 



2 



27rm*ea 3 



(59) 



and the pinning frequency uj p = j^-R c 1 which as 

shown earlier, sets the scale for the conductivity and 
other dynamical quantities, satisfy the following relation 



Rc 
V2a 



(60) 



where, the Larkin length R c ~ R a [£,o I a] 3 ■ Since most 
experiments on the 2D electron gas in Si MOSFETS and 
clean hetero junctions probe dc quantities, they cannot 
permit a simple extraction of the disorder scales. To use 
the available data on dc transport to extract the disorder 
scales, would require a detailed theoretical study of finite 
temperature dc transport, which is clearly beyond the 
scope of this paper and that of the variational method. 
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The only other possibility is through measurements of 
microwave conductivity which presents a straightforward 
way of extracting the pinning frequency. These measure- 
ments in conjunction with SAW measurements which can 
be used to obtain the effective elastic modulus can then 
be combined with theoretical studies to verify the collec- 
tively pinned nature of the crystal. 

To conclude, in this paper, we have studied the the 
physical properties of the Wigner crystal of a low den- 
sity electron gas that is pinned by weak gaussian dis- 
order. We have used the standard combination of the 
gaussian variational method and replicas to calculate a 
host of transport properties. We presented results for the 
evolution of the dynamical conductivity as a function of 
density in both the regimes where the effective particle 
size £ is bigger or smaller than the disorder correlation 
length Tf, the SAW attenuation and the putative power 
radiated by the pinned crystal to its environment. We 
hope that our results provide a stimulus for microwave 
conductivity measurements which have been instrumen- 
tal in establishing the crystalline nature of the 2DEG in 
very strong magnetic fields^. This will then instigate a 
real debate as to whether the insulating phase seen in sys- 
tems which exhibit the 2d MIT is really a Wigner crystal 
or some other phase. 

Clearly, many a question needs to be addressed in order 
to understand the real nature of the crystalline state. As 
mentioned earlier, one first needs to evaluate the elastic 
constants of the crystal as a function of the density. The 



variational wave function used in the Sec. Ill Al to obtain 
the effective particle size can be used to obtain reasonable 
estimates of these constants. Wc expect the bulk mod- 
ulus to not differ much from the result of Ref. |25| since 
this is fixed solely by the coulomb interaction. However, 
the shear modulus is bound to deviate from the classical 
value especially as it approaches melting when the den- 
sity is varied. These values when plugged into the results 
obtained in this paper can lead to a behavior, not too far 
from the melting of the crystal, which deviates strongly 
from that expected for a system with the classical elastic 
constants. Secondly, in the absence of a strong magnetic 
field, the spinful nature of the electrons might lead to in- 
teresting spin physics in the crystalline phase. Interest- 
ing spin liquid behaviors have been seen in studies of the 
multi- exchange spin models in a triangular lattice* and 
it would be pertinent to ask if such exotic spin physics 
arises in the crystal and how does this change when the 
system is pinned by the incipient disorder. Moreover, any 
theoretical study of the magneto-resistance needs to take 
the spin physics into account. Other questions concern 
the behavior of the system close to depinning as also the 
dc transport at finite temperature. Some of these ques- 
tions will be addressed in future work. 
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APPENDIX A: RPA 

A standard way to compute the compressibility of a 
charged system is through the free energy F of the system 
assuming neutrality, and from that to compute dfi/dN. 
Since the system is always supposed to remain neutral 
the compressibility stays finite even for a charged system. 
However, this method is cumbersome since the calcula- 
tion of the free energy is usually more difficult than that 
of a correlation function. In this appendix we show, us- 
ing the RPA approximation, that the formula H23|) indeed 
yields to the standard results for the compressibility of a 
normal metal. 

Using 122fl it is easy to check that the susceptibilities 
Xa/3 = (Papp) are given, in RPA, by 



Xn 

Xl2 



X°Vn x °v 12 

X°Vi2 X°Vi 
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Xl2 



X 




(Al) 



where x° is t ne bare (i.e. for H° only) density-density 
correlation function in one of the systems and the q and uj 
dependence is implicit. It is easy to solve IjAljl to obtain 
the (q dependent) capacitance 



Xii(^,q) - Xi 2 (w,q) 



xV,q) 



l + X°(w,q)(^n(q)-^ia(q)) 

(A2) 



The Fourier transform of the Coulomb potentials are 
given by l|31l) . The true capacitance corresponds to u> = 
and the limit q — > which leads to 



C = 



1 



(2 X Q {q = 0))- 1 +4lTd 



(A3) 



The total capacitance is thus the sum of a geometrical 
one C gcom and one due to the electron gas residing in the 
planes C e \ 



1 



1 



c a 



gcom 



l 



(A4) 



The geometrical one is the standard 1 / (And) result. For a 
simple electron gas, x°(9 = 0) is simply the screening 
length A. One thus recovers that the geometrical dis- 
tance d between the planes is effectively enhanced by the 
electronic screening length A on each plane. 



APPENDIX B: CALCULATION OF THE 
ELASTIC COEFFICIENTS 

The decomposition of the density ( | 10l) allows to re- 
cover quite simply the formulas of Ref. |25| for the elastic 
coefficients. Indeed if one assumes that the interaction 
between the particles of the crystal is 



V(r-r')[p(r)-p }[p(r')-p } (Bl) 



where V(r — r') is the (three dimensional) Coulomb in- 
teraction. The decomposition of the density 1|1U|) allows 
us to rewrite (IB II) as 



H = § J^ Vir - r')[V ■ u(r)][V ■ u(r')] + 

V(r - r ') e 4K -( r - r ') e - lK -( u ( r )- u ( r ')) (B2) 



Pa \ -> 
2 2-f 

K#0 



where we have only kept terms that are not averaged to 
zero due to the translational invariance. The first term in 
(IB2|) gives back the term proportional to q in JSJ . Indeed 
one obtains 



-y 



2 2 

e • Pol 



■in 



u L {q)u L {q) 



(B3) 



This gives the value d = p\jt. 

The second term can be expanded to second order in u 
to give the elastic coefficients. The first order term van- 
ishes because the perfect lattice (u = 0) is an extremum 
of the energy. Up to second order in u, 



Ho 



Po 
2 



V(r 



E 

(u a (r) - u a (r')){u (r) - u (r')) (B4) 



where the summation on the coordinate indices a, (3 = 
(x,y) is implicit. Using 

«a(r) - u a (r') = i u «(q)( eiq r " e4q r ') (B5) 
q 

One can rewrite l|B4|) as 

K#0q,q' Jr > r ' 

« Q (q)«/3(q')(e lq - r - e^ r ')(e lq '' r - e !q '' r ') 

(B6) 
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Using center of mass R = r + r' and relative ro = r — r' 
coordinates can be rewritten as 



i?2 = SE V(r )e^K a K p u a ( q )M-^) 
(2-2cos(q-r )) 

(B7) 

which using the Fourier transform of the potential V(r ) 
denoted U(q) = l/q one obtains 



H * = m EE K a K (2V(K) - V(K + q) - V(K - q)]« a (q)^(-q) (B8) 
q k#o 

I 



which gives the classical elastic constants. 



